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We  analyze  the  laser  cooling  of  a  gas  of  polarizable  particles  by  continuous  dispersive  position  detection  and 
active  feedback.  Microkelvin  temperatures  can  be  attained  inside  an  optical  resonator,  while  in  free  space 
cooling  requires  wavelength-size  beams.  The  maximum  cooling  flux  is  set  by  the  thermal  Doppler  width,  with 
typical  values  between  10^  and  10®  molecules  per  second. 
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A  dissipative  force  on  a  particle  can  arise  from  an  other¬ 
wise  conservative  potential  in  combination  with  a  time  delay. 
For  instance,  in  polarization  gradient  cooling  of  atoms  [1], 
delayed  optical  pumping  between  internal  states  results  in  a 
friction  force.  In  stochastic  cooling  in  particle  accelerators 
[2],  an  electronic  delay  between  a  position  measurement  and 
a  correction  force  is  introduced,  and  in  cavity  cooling  [3,4] 
an  optical  resonator  provides  the  delay.  Electronic  feedback 
cooling  has  been  applied  to  classical  [5]  and  quantum  sys¬ 
tems  [6-10],  and  theoretical  concepts  have  been  extended 
into  the  quantum  domain  [10-13].  For  accelerators,  an  opti¬ 
cal  version  of  stochastic  cooling  has  also  been  proposed  [14]. 

The  feedback  cooling  of  an  atom  strongly  coupled  to  an 
optical  resonator  has  been  analyzed  [15],  and  in  cavity  ex¬ 
periments,  single  atoms  have  been  trapped,  but  not  cooled, 
by  applying  feedback  [16,17].  In  free  space,  various  cooling 
schemes  based  on  optical  measurement  and  feedback 
[18-20]  have  been  proposed,  and  in  an  experiment,  a  sam¬ 
ple’s  center-of-mass  oscillation  has  been  damped  [21]  by 
feedback.  A  method  using  an  optical  resonator  similar  to  the 
one  analyzed  below  [22],  has  been  numerically  studied  by 
Lu  and  Barker  [23].  However,  with  the  exception  of  Refs. 
[10,15,22],  most  treatments  do  not  properly  include  the  pho¬ 
ton  recoil  heating  that  necessarily  accompanies  the  optical 
measurement  and  feedback.  As  discussed  below,  this  can 
lead  to  erroneous  conclusions  regarding  the  feasibility  of  the 
proposed  method  for  a  given  detection  geometry. 

Here  we  analyze  the  laser  cooling  of  a  gas  of  polarizable 
molecules  by  continuous  dispersive  position  measurement 
and  optical  feedback  [22,23].  We  show  that  the  dissipative 
force  is  proportional  to  the  radiation  pressure  associated  with 
Rayleigh  scattering  into  the  detector,  while  its  velocity  de¬ 
pendence  is  determined  by  the  frequency  dependence  of  the 
loop  gain.  The  maximum  collective  cooling  rate  for  a 
sample,  limited  by  collective  heating,  is  given  by  the  thermal 
Doppler  effect.  This  results  in  a  typical  cooling  flux  of 
10^  molecules/s  at  room  temperature,  and  10^  molecules/s 
at  a  sample  temperature  of  100  ptK.  Using  a  resonator  for 


signal  enhancement,  the  method  can  be  applied  to  any  gas  of 
particles  that  scatter  light  at  a  sufficiently  large  rate. 

The  dissipative  force  arises  from  the  application  of  a 
time-dependent  conservative  optical  potential  in  response  to 
an  optical  measurement  of  the  particle’s  motion  [22,23]. 
Consider  a  single  molecule  or  atom  of  mass  m  moving  as 
x„,=vt  at  constant  velocity  u  in  a  weak  periodic  potential 
F(x,f)  =  {/(?)cos  {2kx)  whose  depth  U{t)<ymv^  depends  on 
the  molecule’s  motion.  Since  the  time  variation  of  the  force 
on  the  molecule  f{t)  =  2kU{t)sm{2kvt)  arises  both  from  the 
spatial  variation  cos  (2kx)  of  the  potential,  and  from  the  de¬ 
pendence  of  the  depth  U{t)  on  the  molecule’s  trajectory,  a 
frequency  component  of  t/(?)  in  phase  with  the  molecular- 
motion-induced  force  variation  sin(2A:i;f)  produces  a  nonzero 
average  force  that  can  heat  or  cool  the  molecule.  A  standing 
light  wave  can  be  used  to  both  monitor  the  position  Xm(f)  of 
the  polarizable  particle  by  means  of  its  index  of  refraction 
[3,4,17],  and  to  generate  a  periodic  potential  of  adjustable 
depth  U{t)  via  the  light  shift. 

For  a  sample  instead  of  a  single  particle,  and  linear  feed¬ 
back,  the  self-generated  dissipative  force  of  each  molecule 
remains  unchanged,  while  the  fluctuations  in  U{t)  due  to 
other  molecules  constitute  a  source  of  heating.  However,  the 
heating  is  quadratic  in  loop  gain,  whereas  the  cooling  varies 
linearly  with  gain  [2].  Therefore  the  collective  heating  limits 
only  the  cooling  rate,  but  not  the  final  temperature  attainable 
at  low  loop  gain. 

The  optical  signal  arising  from  the  particle’s  motion,  and 
consequently  the  cooling  force,  can  be  increased  by  means  of 
an  optical  resonator  (Fig.  1).  Light  of  angular  frequency  ck 
incident  onto  a  standing-wave  resonator  of  finesse  F,  waist 
size  w,  and  field  decay  rate  constant  y^,  produces  an  electric 
field  of  amplitude  2E^  at  an  antinode  (x=0)  inside  the  reso¬ 
nator.  A  molecule  with  complex  polarizability  a  moving  on 
the  resonator  axis  coherently  scatters  photons  into  free  space 
at  an  average  rate  rs;.=^^|Re(a)£'^p/(67reo^),  and  experi¬ 
ences  an  optical  potential  Uq  cos  {2kx)  of  depth  Uq 
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FIG.  I.  (Color  online)  Schematic  of  continuous-feedback  cool¬ 
ing.  The  moving  molecule  periodically  detunes  the  standing-wave 
resonator  by  S„,{t).  The  observed  fractional  change  in  cavity  trans¬ 
mission  E(f)  acts  as  an  error  signal  that  is  applied  to  change  the 
input  power  P,  by  a  fraction  -Ef|,(r).  The  feedback  open-loop  gain 
H(s)  is  set  by  an  external  electronic  circuit. 


=  -|£’cpRe(a')/2.  Conversely,  the  forward-scattered  field,  or 
equivalently,  the  molecule’s  index  of  refraction,  results  in  a 
molecular-position-dependent  detuning  Si„  =  ^y^  cos{2kx,„)  of 
the  resonator  from  its  average  resonance  frequency  [3,4]. 
Here  (,  =  iir]T Uq  is  a  dimensionless  parameter  characteriz¬ 
ing  the  molecule-cavity  coupling,  and  the  single-atom  coop- 
erativity  2ri=\2F ! {Trk^vP')  can  be  interpreted  as  the  fraction 
of  photons  scattered  into  the  resonant  cavity. 

For  incident  light  detuned  by  an  amount  Sj=ck-cOc  rela¬ 
tive  to  the  resonator,  the  time-varying  detuning  of  the 
resonator  by  the  moving  molecule  changes  the  resonator 
transmission  by  a  fraction  e(f).  The  transmitted  power  is 
measured  and  used  as  an  error  signal  in  a  feedback  loop  to 
adjust  the  incident  power  by  a  fraction  -efb(t)  (Fig.  1).  In 
order  to  ignore  cavity-induced  forces  arising  from  the  finite 
response  time  of  the  resonator  [3,4],  we  take  its  linewidth 
2  to  be  much  larger  than  both  the  feedback  bandwidth  and 
the  Doppler  frequency  2kv,  such  that  the  intracavity  power 
P^(f)  adjusts  instantaneously  to  a  value  determined  by  the 
total  light-resonator  detuning  S,{t)  =  Si- .  If  both  the 
molecule  and  the  feedback  loop  produce  only  small  changes 
(|5„,|/yc>l®fbl  ^  then  the  fractional  deviation  e(f) 
=  Pc{t)l Po-  \  of  the  intracavity  power  Pjj)  from  its  unper¬ 
turbed  value  Pq  (for  <5m  =  6fb=0)  is 

e(r)  =  r7;'4„(f) -eft,(t).  (1) 


The  moving  molecule  modulates  the  intracavity  power  by  an 
amount  rS„J  proportional  to  the  normalized  resonator 
slope  r=2Siy^l {'^^  +  ^),  while  the  feedback  loop  adjusts  the 
incident  power  by  a  fraction  -efi,(f).  If  the  molecule’s  kinetic 
energy  far  exceeds  the  light  shift  Uq,  then  to  lowest  order 
is  determined  by  the  molecule’s  unperturbed  motion 
x,„=vt.  We  introduce  the  open-loop  feedback  gain  H{s)  for 
the  Laplace  transformed  quantities  ef|3(j),  8(5)  via  ef^,=Hs. 
Then  +  with  real  and  imaginary  parts 

Hi  and  H2,  respectively,  is  the  complex  gain  in  the  frequency 
domain.  If  the  loop  is  stable,  the  steady-state  solution  of  Eq. 
(1)  in  the  time  domain  is  given  by 


e(r)  = 


[1  +  Hi{2kv)]cos{2kvi)  +  H2i2ku)sm{2kvt) 
\1+Hi2ikv)\^ 


(2) 


Hi{2kv)  and  H2{2kv)  are  the  open-loop  gain  in  phase  and  in 
quadrature  with  the  molecular-motion-induced  intensity 
modulation  S„^{t)  =  ^y^cos{2kvt),  respectively.  They  deter¬ 


mine  the  closed-loop  signal  e(f),  and  thereby  the  time  varia¬ 
tion  of  the  optical-potential  depth  {/(r)  =  I7o[l+e(t)].  In  the 

limit  Uq<2'™  ^  work  is  done  on  the  molecule  at  a  rate  W 
=  e{i)fjj)v  to  lowest  order,  where /„=2A:I7o  sin(2/cut)  is  the 
unperturbed  force.  The  component  sin(2A:ur)  of  e(f)  in  phase 
with  /„  produces  a  dissipative  force  with  time  average  / 
=  {W)lv  given  by 


f{v)=fikr]T^^ 


rH2{2kv) 

\+H{2ikv)\^' 


(3) 


This  expression,  valid  for  arbitrary  laser  detuning  from  mo¬ 
lecular  resonances  below  saturation,  shows  that  the  friction 
force  /  depends  on  molecular  parameters  exclusively  through 
the  Rayleigh  scattering  rate  Fs^,.  /  is  proportional  to  the  rate 
of  momentum  transfer  2tikr]T^^  due  to  backward  scattering 
into  the  resonator  at  rate  multiplied  by  a  dimensionless 
function  of  the  molecule  velocity.  In  particular,  the  sign  and 
velocity  dependence  of  /  are  determined  by  the  frequency- 
dependent  loop  gain  H{ico).  The  dissipative  force  /  is  maxi¬ 
mum  for  a  resonator-light  detuning  <5,  =  ±7^  that  gives  the 
largest  slope  r=±l.  In  the  following  we  assume  4  =  “7c  ('' 
=-l),  such  that  H{s)>0  corresponds  to  negative  feedback. 

The  cooling  force  /  is  proportional  to  the  quadrature  com¬ 
ponent  sm{2kvt)  of  the  intracavity  light  modulation  in  closed 
loop,  given  by  //2/|l-t//p,  evaluated  at  the  Doppler  fre¬ 
quency  2kv.  In  particular,  for  very  small  or  very  large  open- 
loop  gain  (|//p^l  or  |//p>l),  this  relevant  quadrature  of 
the  intensity  variation,  and  hence  /,  will  be  small.  The 
velocity-dependent  term  takes  on  its  maximum  value  1/(2 
+  2Hi)  when  the  quadrature  gain  H2  and  the  in-phase  gain  Hi 
are  related  by  I//2I  =  |  1  When  |1h-//i|^  1,  the  feedback 
loop  regeneratively  amplifies  the  intensity  variation  caused 
by  the  moving  molecule.  However,  the  heating  of  the  mol¬ 
ecule  due  to  noise  amplification  is  then  also  increased. 

The  simplest  stable  cooling  loop  is  a  differentiator  [17] 
H^{io))  =  ico/{2ku),  with  unity  gain  frequency  2ku.  (The  gain 
can  be  rolled  off  outside  the  velocity  range  of  interest  to  the 
cooling.)  We  assume  for  simplicity  that  the  laser  is  far  de¬ 
tuned  from  molecular  resonances  compared  to  the  Doppler 
shift  2kv,  such  that  F^c  is  independent  of  v.  The 
differentiator-induced  cooling  force  is  then  given  by 


fdiv)  =  -  hkrjT,^—^ - 5 .  (4) 

u  +  V 

The  maximum  force  fikr]T^J2,  attained  for  the  unity-gain 
velocity  u  =m,  is  the  same  as  in  conventional  Doppler  cooling 
at  the  photon  scattering  rate  77Fs;./2,  but  the  differentiator 
force,  falling  off  as  1/u  rather  than  1/u^,  provides  a  signifi¬ 
cantly  larger  velocity  capture  range.  Figure  2  shows  that 
higher-order  loops  can  substantially  extend  the  velocity  cap¬ 
ture  range  while  maintaining  the  low-velocity  friction  coef¬ 
ficient  aQ=dfldv  that  determines  the  final  temperature. 

Technical  or  quantum  noise  8(r)  in  the  light-induced  po¬ 
tential  [/(t)  =  I7o[l+e(t)]  will  heat  a  particle  moving  at  ve¬ 
locity  V  in  proportion  to  the  spectral  noise  density  at  the 
modulation  frequency  2kv  of  the  unperturbed  force.  We  de¬ 
fine  the  single-sided  fractional  spectral  noise  density 
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FIG.  2.  (Color  online)  Cooling  force  versus  molecular  velocity  v 
for  conventional  Doppler  cooling  (a),  and  for  feedback  cooling  with 
loop  gain  H{s)=s  (differentiator,  b),  //(^)  =  i(l+i/10)/(l-l-8i/10) 
(c),  and  //(i)  =  i(l-Hi/10)(l-Hi/100)/[(l-t8i/10)(l-l-i/20)]  (d). 
Here  s  =  iv/u,  where  u  sets  the  unity  gain  frequency  for  the  differ¬ 
entiator  loop. 

=  (2/7r)/gC?t(8(0)e(f))cos(wr),  normalized  such  that 
/Qt/w5'c((u)=(e^(0))  is  the  mean-square  fractional  intensity 
noise,  and  calculate  the  heating  power  due  to  the  force  fluc¬ 
tuations  [24], 

Trk^ 

- UlS,{2kv).  (5) 

m 


photoelectron  detection.  In  addition  to  the  cavity  heating 
VLsn,  scattering  into  free  space  heats  the  molecule  at  a  rate 

IVfs  =  2£’rrs;,  for  isotropic  scattering  [24]. 

For  the  differentiator  loop  Hj^ici))  =  io)l{2ku),  the  choice 
u  =  {q~^ -  \+2rf^)fikl m  for  the  unity-gain  velocity  u  mini¬ 
mizes  the  equilibrium  temperature  Td,  yielding 

^Br,  =  4£,(77-'  +  l)(277-'  +  g-i-l),  (7) 

where  kg  is  Boltzmann’s  constant.  The  temperature  scales 
with  the  recoil  energy  E,.,  but  strongly  depends  on  the  coop- 
erativity  77,  i.e.,  the  molecule-cavity  coupling.  In  the  bad- 
cavity  limit  r]<q<  \  the  final  temperature  is  independent  of 
detector  quantum  efficiency  q,  and  given  by  Tj=%rf^E^lkg. 
For  a  typical  bad-cavity  setup  with  77=0.05  for  Cs  atoms 
[25],  this  corresponds  to  a  temperature  of  320  /xK,  while  the 
good-cavity  limit  77^  1  [16,17]  with  <7  =  0.2  would  corre¬ 
spond  to  Tj=\.6  /U,K. 

For  a  thermal  sample  consisting  of  N+ 1  molecules  at 
temperature  T=mv]jkg,  any  chosen  probe  molecule  at  ve¬ 
locity  V  will  be  heated  by  the  intracavity  intensity  noise  in¬ 
duced  by  the  other  N  randomly  moving  molecules.  Substitut¬ 
ing  the  corresponding  closed-loop  spectral  noise  density  into 
Eq.  (5)  yields  the  collective  heating  rate  W[^  as  follows; 


This  general  formula  describes  the  heating  of  a  molecule 
moving  at  velocity  v  due  to  classical  intracavity  intensity 
noise  with  spectral  noise  density  S^io)).  For  quantum  noise, 
Eq.  (5)  needs  to  be  doubled  to  take  into  account  atomic- 
dipole  fluctuations  [24].  Neglecting  technical  noise,  the  pho¬ 
ton  shot  noise  of  the  incident  light  with  power  P,  corre¬ 
sponds  to  a  fractional  power  spectral  density  given  by  5]’“ 
=  hck/{PiTr).  In  closed  loop,  the  intracavity  noise  density  is 
5P’^‘'(a7)  =  5[’“/|l  H-//(fw)p  for  noise  frequencies  (o<  In  ad¬ 
dition,  uncorrelated  shot  noise  in  the  detector  photocurrent  at 
limited  quantum  efficiency  q<l  arising  from  the  undetected 
photons  will  cause  the  feedback  loop  to  generate  (classical) 
intracavity  intensity  noise.  This  results  in  a  fractional  noise 
density  5[™=5'[’“(q'"^- l)|//p/|  1  in  closed  loop.  The  to¬ 
tal  shot  noise  expressed  in  terms  of  the  intra¬ 

cavity  power  P^=P^F/{2tt)  for  the  detuned  cavity  with  |4| 
=  y^,  and  inserted  into  Eq.  (5)  then  yields  a  closed-loop  heat¬ 
ing  due  to  shot  noise 


= 'E,77r,, 


2+\H{2ikv)Wq-^-l) 
|1  +Hi2ikv)\^ 


(6) 


where  E^=hV !  {2m)  is  the  recoil  energy.  We  see  that  Wsn 
can  be  expressed  as  the  recoil  heating  rate  2Ei.r]T^^,  multi¬ 
plied  by  a  dimensionless  function  of  the  loop  gain  El  and 
detector  quantum  efficiency  q.  In  the  absence  of  feedback 
{H=Q),  Eq.  (6)  reproduces  the  well-known  recoil  heating 
[24]  due  to  scattering  into  the  resonator,  where  each  back- 
scattering  event,  occuring  at  a  rate  {rf/ 2)rs^  for  the  detuned 
cavity,  heats  the  molecule  by  an  amount  4£^.  Eor  unity  quan¬ 
tum  efficiency  q=l,  the  feedback  loop  reduces  (increases) 
this  heating  by  suppressing  (enhancing)  the  intracavity  fluc¬ 
tuations  for  |1h-//|>1  (|1h-//|<1).  For  q<l,  the  feedback 
loop  also  generates  intensity  noise  originating  from  random 


1+Hi2ikv)\'^  2kuth 


(8) 


As  is  typical  of  stochastic  cooling,  the  heating  overwhelms 
the  cooling  at  too  large  sample  size  N  or  too  large  photon 
scattering  rate  Fs^,.  The  collective  heating  rate  is  quadratic 
rather  than  linear  in  Fs^,,  and  can  be  reduced  below  the 
single-particle  heating  at  the  expense  of  cooling  speed.  The 
cooling  rate,  however,  is  limited  by  the  collective  heating. 
For  a  differentiator  loop,  the  net  cooling  power  averaged 
over  the  thermal  sample,  {W^i)={f^v  +  Wf^)  is  maximized  for 
u~v±.  The  cooling  rate  constant  y^=-2{W^ ! {kgT)  can  then 
be  written  in  the  form 


y,«^(2f_-fL),  (9) 

oJy 

where  T^^^=2Nr)T^Jil{kgT)  is  the  total  scattering  rate  into 
the  cavity  normalized  to  the  sample  temperature.  The  cooling 
rate  is  maximized  for  r;.av=l,  yielding  a  rate  constant  y^ 
=  kvi'al  {bN).  We  see  that  the  thermal  Doppler  broadening  kv^^ 
takes  the  role  of  the  stochastic-cooling  bandwidth  [2],  and 
that  the  cooling  of  smaller  subsamples  can  proceed  faster 
[2,18].  For  a  sample  of  A=10^  magnetically  trapped  CaH 
molecules  at  0.4  K  [26]  cooled  with  light  of  wavelength 
760  nm,  the  optimum  cooling  rate  is  7^=0.!  s"',  attained  at 
a  photon  scattering  rate  77Fs;,=3  X  10^  s”'.  (Here  we  are  as¬ 
suming  that  the  different  degrees  of  freedom  are  mixed  [2], 
e.g.,  by  the  trapping  potential  [18].)  An  ensemble  of 
10^  CaH  molecules  at  room  temperature  can  be  cooled  at  a 
rate  7^=340  s”'  for  a  photon  scattering  rate  of  77Fsj,=3 
X  10^  s“'.  Figure  3  shows  the  differentiator  cooling  flux  yJSl 
for  different  sample  sizes,  ranging  between  10^  and 
10®  molecules /s  for  temperatures  between  100  mK  and 
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FIG.  3.  (Color  online)  Differentiator  cooling  flux  (cooling 
power  divided  by  kgT)  versus  temperature  T  for  various  free-space 
scattering  rates  and  sample  sizes  N,  calculated  for  CaH  with  X 
=  760  nm  inside  a  bad  cavity  (77=0.05).  The  parameters  are  (a)  N 
=  10®,  rsc=10^s“'  (thin  solid  line),  (b)  Al=10^,  rsc=10^s“* 
(dashed),  (c)  A^=10®,  rs,,=  10^s“*  (dotted),  and  (d)  Af=10*, 

=  10  s“*  (dash-dotted).  The  thick  solid  line  is  the  optimum  rate  with 
Fs,,  adjusted  to  maintain  F|,2^,=  l. 

1000  K.  This  is  somewhat  lower  than,  hut  comparable  to,  the 
atomic  fluxes  attained  with  standard  laser  cooling  methods 
[27]. 

The  prediction  of  heating  for  too  large  cooling  speed 

Fcav  >  2  by  Eq.  (9)  is  not  necessarily  correct.  The  reason  is 
that  self-organization  of  the  atom-light  system  [25]  may  oc- 
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cur.  The  corresponding  change  from  single-atom  Rayleigh 
scattering  to  much  stronger  Bragg  scattering  by  the  sample 
could  then  lead  to  new  collective  effects  not  captured  by  the 
above  model. 

The  role  of  the  resonator  is  to  enhance  the  photon  scatter¬ 
ing  rate  77rsj,  into  the  detector,  i.e.,  the  signal-to-noise  ratio 
[19].  In  the  absence  of  the  resonator,  rj  is  replaced  by  the 
detection  solid  angle  Afl  =  3/(A:^H'^),  preventing  cooling  alto¬ 
gether  except  for  the  smallest  beam  size.  For  instance,  the 
cooling  limit  for  Rb  atoms  using  a  beam  as  small  as  w 
=  10  /zm,  according  to  Eq.  (7)  corresponds  to  7^=370  K,  i.e., 
no  observable  cooling.  This  explains  why  stochastic  cooling 
has  not  been  observed  in  free-space  experiments  using  mac¬ 
roscopic  beam  size  for  detection  [18,21]. 

In  conclusion,  we  have  derived  simple  analytic  expres¬ 
sions  for  the  cooling  and  heating  of  a  gas  of  linearly  polar¬ 
izable  particles  interacting  with  a  laser  beam  whose  intensity 
is  adjusted  in  response  to  the  particles’  motion.  The  velocity 
capture  range  can  be  tailored  via  the  frequency-dependent 
external-loop  gain,  while  the  thermal  Doppler  width  takes 
the  role  of  the  stochastic-cooling  bandwidth,  determining  the 
maximum  cold-molecule  flux. 

This  work  was  supported  by  the  NSF,  DARPA,  and  the 
Sloan  Foundation.  J.S.  acknowledges  NDSEG  and  NSF. 
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